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Abstract 

We analyzed the solar, terrestrial and atmospheric neutrinos experiments us- 
ing the three- flavor neutrino framework and got the allowed regions for param- 
eters (Amf 2) sin 2 2^12, Am^, #13, #23)- In solar neutrino experiments, we got 
the large angle solution (Amf 2 , sin 2 26 12 ) = (4x lCT 6 -7x lCT 5 eV 2 , 0.6-0.9) 
and small angle solution (3 x 10~ 6 - 1.2 x 10~ 5 eV 2 , 0.003 - 0.01) for 
#13 = 0° — 20°. From the terrestrial and atmospheric neutrino experiments, 
we got the allowed regions (6 13 < 4°, 24° < 6» 23 < 26°) for Am^ 3 = 2eV 2 , 
(6»i3 < 4°, 24° < 6*23 < 45°) for Am| 3 = 0.2eV 2 and (6 13 < 14°, <9 23 ~ 40°) for 
Am^ = 0.02eV 2 . It seems that the large angle solution for (Am 2 2 , sin 2 26*12) 
is favored than the small angle slution from the analysis of zenith angle de- 
pendence in atmospheric neutrino sub-GeV experiment. 



I. INTRODUCTION 

The problem of neutrino masses and oscillations is one of the most interesting issues 
to study physics beyond the standard model (SM) [IJ. In many experiments which are 
under way, indications in favor of neutrino masses and oscillations have been obtained. In 
these experiments, the solar neutrino experiments |@-45| measure the event rates significantly 
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lower than the ones predicted by the standard solar model, and the atmospheric neutrino 
experiments [0-0] observe an anomaly referred to as the unexpected difference between the 
measured and predicted /i-like and e-like neutrinos. Another indication in favor of non-zero 
neutrino masses is in the cosmological analysis by dark matter [[J. 

On the other hand, terrestrial neutrino experiments searching for the neutrino masses 
and oscillations are under way. The E531 [TO], CHORUS and NOMAD [I I experiments 
using the beam from accelerator search for v T appearance in a and E776 | I2]] , KARMEN 
|T3| and LSND |TJ[ experiments using the accelerator beams are searching for — > v e 
and — > v e oscillations. The experiments using nuclear power reactor |D| search for 
the disappearance of u e , in which v e — > i>x {X — //, r) transitions are expected. These 
experiments do not observe significantly large neutrino transitions. 

In this paper, we analyze the data of current solar, terrestrial and atmospheric exper- 
iments in a framework where the three neutrinos have masses and mix each other, and 
search the allowed regions of parameters (Am^, sin 2 2#i 2 , Am| 3 , 6*13, #23) characteriz- 
ing the masses and mixing of three-flavor neutrinos. Although there are many analyses 
which study the solar, atmospheric and terrestrial neutrino problems in three-flavor neutri- 
nos framework ||i6|--|T8|j , we study more thoroughly these problems including the analysis of 
zenith angle dependence in recent Super Kamiokande atmospheric experiment p!9|j20|| . After 
the analyses of experiments for neutrino oscillation, we will present a matrix of neutrino 
mixing ascribed from the allowed regions of parameters obtained. 



II. NEUTRINO OSCILLATION 

Weak currents for the interactions producing and absorbing neutrinos are described as 
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h = 2 J2 iLalftU^pVLp, (1) 
a, 0=1 

where l a (h = e, I2 = //, h = t) represents the lepton flavor, vp the neutrino mass eigenstate 
and U is the lepton mixing matrix. U is the unitary matrix corresponding to the CKM matrix 
Vq KM for quarks defined by 

U = UiUl (2) 

where the unitary matrices Ui and U v transform mass matrices Mi for charged leptons and 
My for neutrinos to diagonal mass matrices as 

UiM l U^ 1 = diag[m e , m M , m T ], 
U u M v Uy 1 = diag [m 1} m 2 ,m 3 }. 

We present the unitary matrix neglecting the CP violation phases as 
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JJ — g«623A7gi6l3A5gifll2A2 



'12^13 



'12 c 13 



'13 



'12^23 u 12°23' : '13 u 12 u 23 °-\2°23 c, 13 °2^1'i 
12^23 ' b 12 c 23 i, 13 c 23 c 13 



'12^23 ^12^23°13 

c& = cos0£, s 



(4) 



sinC 



where Aj's are Gell-Mann matrices. 

The probabilities for transitions v\ a — > v\ are written as 



P(yi a -> %)= I < %WK(o) > | 2 = <f w . + Pi 



12 



23 



S23 + V 



31 



5, 



- -n«-«^3i, 



- -2* w ,(l - 2f/ 2 3 ) + 2(0^0* + 0* 2 L/* - C/^C/f 



p 23 

31 



-25 lah {\ - 2Uf al ) + 2(-Ul l U? 0l + Ul 2 Uf B2 + U&Jt 
-25 lalp (l - 2UfJ + 2(0* ^ - C/f a2 C/^ 2 + IJj;., 



P 



j2 r 7-2 

t2 tt2 



where Sy is the term representing the neutrino oscillation; 

9 Am. 2 , 
% = sin 2 1.27— ^L, 



(5) 



(6) 



in which Am^ = |w| — m||, and L are measured in units eV 2 , MeV and m, respectively. 
From the unitarity of U, we get relations 

P% a ^e + + V% a ^ T = 0, z, j = 1, 2, 3, / Q = e, ji, r. (7) 

The values of neutrino masses are not known precisely, but we know that if one identifies 
the dark matter of universe (or at least its hot dark matter component) with neutrino matter 
one has [01 



m\ + ni2 + 1713 ~ several eV. 



(8) 



We do not use this value strictly in present analysis, but we consider that the sum of neutrino 
masses is not so small. In two-flavor neutrinos analyses in which one mass parameter Am 2 
appears for solar neutrino experiments, the result that Am 2 is 10~ 4 — 10 _5 eV 2 or ~ 10~ 10 eV 2 
is obtained [ 21 1 . For atmospheric experiments, Am 2 is obtained as 10" 1 -10- 2 eV 2 lEl. Then 
it seems that two neutrinos masses in three neutrinos are very close and another one is rather 
far away from them. Then we assume that three neutrino masses have such a mass hierarchy 
as 



mi ~ <C 1713. 



(9) 



In the the mass hierarchy Eq. Am 2 2 <C Am 2 3 ~ Am 2 3 , the expression Eq. (5) for the 
transition probabilities P{yi a — > Pu) are rewritten as 
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P(u e -> u e ) = 1 - 2(1 - 2f/ 2 3 - C/ 4 - f/ 4 2 + t/ e 4 3 )5 12 - 4f/ 2 3 (l - 2 3 )S 23 , 
P(i/„ -> = 1 - 2(1 - 2f/ 2 3 - Ift - 4 2 + 4 3 )S 12 - 4C^g(l - C/ 2 3 )S 2 3, 
P(i/t _> Vt ) = i _ 2(1 - 2f/ 2 3 - £ft - f/ 4 2 + f/ 4 3 )S 12 - 4I&(1 - f/ 2 3 )S 23 , 
P(i/„ - i/.) = P(z/ e - = 2(E^E^ + [/ 2 2 C/ 2 2 - U%U^)S 12 + 4C/ e 2 3 C/ 2 J ^ 3 , 
P(i/ T - v e ) = P{v e -> i/ T ) = 2(0* L/* + 2 2 e 2 2 - r 2 3 e 2 3 )S 12 + 40 2 3 2 3 S 23 , 
P(z/ T - „„) = P(u, - „ T ) = 2(C/ 2 1 C/ 2 1 + 0* 2 [/ 2 2 - U%J%)Su + ^ 3 U 2 r3 S 23 . 

(10) 



III. NUMERICAL ANALYSES OF NEUTRINO OSCILLATION IN 
THREE-FLAVOR NEUTRINOS 

A. Solar neutrinos 

We first analyze the solar neutrino experiments. Considering the matter effects (MSW 
effect [p2" D in three-flavor neutrinos, the transition probability for v e — > v e is expressed as 
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P 3 ^ sw (Am 2 2 , 12 , 13 , P) = cos 4 ^ 13 P 2 ^ sw (Am 2 2 , 12 , i3 , P) + sin 4 0i 3 C (11a) 

P 2l ! s w (Am 2 2 , 6 12 , 13 , P) = ^ + Q - 6(A cos 2 13 - Am 2 12 cos 20 12 )P C (0 12 , P) 

x cos 20 i2 cos 202, (Hb) 
P c (0i 2 , P) = exp(-^ 7 (0 12 ,P)), (11c) 

fa t?\ Am 2 2 sin 2 20i 2 

7(0i 2 , P = — — , lid 

2P cos 20i2 1 dN e /N e dx \ 

where G is the theta function, A = 2\/2GpN e E {Gp is Fermi constant, N e number of electron 
per cm 3 at the production point of neutrinos in the sun and P the energy of neutrino), P c 
the Landau-Zener-Stueckerberg crossing probability, 7 is the adiabaticity parameter (| • • • | 
represents the value at the production point of neutrinos) and 0^ the mixing angle at 
the production point. P 2 ^ sw is the transition probability with a replacement of N e — ■> 
N e cos 20i 3 in two-flavor neutrinos transition probability. This expression is obtained from 
an approximation; 

/ 9 ./x „ 1 ( * (mi -A/2) sin 2 20i 3 \ , . 

G[A - m 2 - A/2 cos20 13 exp -- , 3 Ar ' — - < 1, 12 

L v 3 11 J V 4 I dN e /N e dx 1 P cos 20 i3 / v ; 

where A = ((m\+mty — (m\—m\) cos20 12 ). This approximation is reasonable for present as- 
sumption of mass hierarchy Eq. (9) because of A -C m|— A/2 and m|— A/2 ^> \dN e /N e dx\oE . 



4 



The ratios of the detected e neutrino fluxes to the expected e neutrino fluxes deduced 
from the standard solar model (SSM) |23[ are expressed by using the transition probability 



pMSW „„ 

= it::p^ sw (E)f(E)dE 

where /(-E) is the product of the spectral function of neutrino flux and detector sensitivity. 
The neutrino flux is the sum of many neutrino fluxes produced by the various nuclear fusion 
reaction at the center of the sun. The detector sensitivity also depends to the neutrino 
energy. Numerical results for neutrino flux and detector sensitivity are analyzed precisely 
in Ref. f23]| . For f(E) used in this paper, see Appendix A. 



In our analysis, we use the following three experimental data for R: 
Ga experiment : 

R = 0.534 ±0.087, (14a) 

CI experiment [ffl: 

R = 0.274 ± 0.046, (14b) 

water Cherenkov experiment |J: 

R = 0.437 ±0.092. (14c) 

In Ga experiment, we combined the SAGE || and GALLEX || data. First we estimate 
the numerical values R for # 13 = which corresponds to two-flavor neutrinos case using the 
Eq. (13), and show the contours of R on sin 2 29u — &-rn\ 2 plane in Fig. |l} Each contour 
denoted as Ga, CI and Kam corresponds to the upper and lower values of R in Eq. (14) 
for Ga, CI and Kamiokande's water Cherenkov experiment. There are two solutions which 
are denoted as common areas enclosed by each two contours of Ga, CI and Kam. These 
are called as large solution and small solution. This result is similar to the one obtained in 
various analyses |24[ . 



Next we estimate the numerical values R for the case of #i 3 7^ 0. In Fig. ^|, we show the 
allowed regions of the combined Ga, CI and Kam experiments using the x- sc i uare > where 
the solid thin, solid thick and dotted lines define the regions allowed at 99%(x 2 = 9.2), 
95%(x 2 = 6.0) and 90%(x 2 = 4.5) C.L., respectively. Fig. |(a) - Fig. |(h) show the allowed 
regions for #13 = 0° — 50°. Fig. |2|(a) shows the #13 = case, thus this shows the two-flavor 



neutrinos' solution. These results are similar to the ones obtained by Ref. [18]. Numerically, 
we show the allowed regions in 95% C.L.; 



for 13 = 0° - 20° 

?„ = 4 x 10~ 6 - 7 x 10~ 5 eV 2 

> large angle solution (15a) 



^'"•12 T: /\ -L W 

sin 2 2^i2 = 0.6-0.9 
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Am 2 12 = 3 x 1(T 6 - 1.2 x l(T 5 eV 2 

sin 2 2^i2 = 0.003 - 0.01 

for 9j3 = 25° - 40° 



> small angle solution 



(15b) 



Am 2 2 = 2 x 10~ 6 - 3 x 10~ 5 eV 2 

sin 2 2# 12 = 0.001 - 0.01 
for #i3 = 45° - 50° 

2 x 10~ 6 - 3 x 10~ 5 eV 2 



> small angle solution 



(15c) 



Am 2 2 



sin 2 26 



(156) 
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0.001-0.7 



The characteristic feature of three-flavor neutrinos' solution is as follows; increasing 8 13 from 
0° to 25°, the small mixing solution and the large mixing solution merge in a single solution, 
further increasing #13 to 50°, the allowed region becomes broader and next shrinks and lastly 
disappears. 



B. Terrestrial neutrinos 



In terrestrial experiments, there are two types of experiment: short baseline and long 
baseline experiment. In present study, we analyze the short baseline experiment. In the 
short baseline experiments, there are E531 flQ|| , CHORUS and NOMAD |[l| accelerator 
experiments searching for v T appearance in v^. We use the data of E531, CHORUS and 
NOMAD experiments; 



P(Vu 



i/ T ) < 2 x 10" 3 (90% C.L.), 
L/E ~ 0.02. 



(16) 



For the experiments searching for v„ — > v e and zAj — > v e oscillations, there are E776 [12 



KARMEN |T3[ and LSND accelerator experiments 



P(u„ 



P(u u 



P{Vn 



u e ) < 3 x 10~ 3 (90% C.L.), E776 
L = 1km, E ~ lGeV, 
v e ) < 3.1 x 10" 3 (90% C.L.), KARMEN 

L = 17.5m, E < 50MeV, 

3 



Ve. 

L 



3.4+ 2 ;° ±0.7 x 10" 3 , LSND 



(17a) 
(17b) 
(17c) 



30m, 



E ~ 36 - 60MeV. 



Furthermore, we analyse the experiments using nuclear power reactor [y| searching for the 
disappearance of v e , 
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1 - P(v e -> u e ) < 1(T 2 (90% C.L.), (18) 
L = 15, 40, 95m, E ~ 1 - 6MeV. 

In the short baseline experiments, the neutrino propagation length L is about L ~ 1km, 

then = sin 2 1.27 — — — L is very small because Amf 2 ~ 10~ 5 — 10 _4 eV 2 suggested from 

E 

last solar neutrino analyses. Then, the £23 term is dominant in the transition probability 
Eq. (10). Seeing that the mixing parts proportional to S23 are 4£7" 2 3 (1— ?7 2 3 ), 4£/ 2 3 (l — £/ 2 3 ) and 
4[/ 2 3 (l — t/ 2 3 ), the transition probabilities in present short baseline terrestrial experiments 
seem to be insensitive to the parameters; Am 2 2 and #i 2 . 

We show the contour plots of the allowed regions on (tan 2 # 13 , tan 2 6*23) plane determined 
by the probability P expressed in Eq. (10) and above experimental data Eqs. (16), (17) and 
(18) in Fig. 0. Allowed regions are corners, left and right hand sides surrounded by curves. 
Curves represent the boundary of 90 % CL. of P. We fixed the parameters Am 2 2 and 
9 12 as Am 2 2 = 10~ 5 eV 2 and sin 2 29 12 = 0.8, and the parameter Am 2 3 to be various values 
from 0.02eV 2 to 20eV 2 . Although we fix the parameters Amf 2 and #12 as Am 2 2 = 10~ 5 eV 2 
and sin 2 29\ 2 = 0.01, the results are not changed. Dotted lines show the allowed regions 
restricted by the LSND data. From these results, we see that the numerical allowed regions 
on (#i3, # 23 ) without the LSND data are as follows; 

for Am 2 3 = 20eV 2 

(<4°, <2°), (<2°, >88°), (>86°-88°, any), (19a) 
for Am 2 3 = 2 eV 2 

(<4°, <26°), (<2°, >65°), (>86°-88°, any), (19b) 
for Am 2 3 = 0.2 eV 2 

(<4°, any), (> 86°, any), (19c) 
for Am 2 3 = 0.02 eV 2 

(< 12°, any), (> 78°, any) (19d) 
for Am 2 3 < 0.005eV 2 

all regions are allowed. (19e) 

If we combine the LSND data with the above analyses, the allowed region disappears 
lower than 0.2 eV 2 of the Am 23 value. Furthermore, combining the solar neutrino solutions 
with this terrestrial ones, the allowed regions larger than 50° of 9 13 on terrestrial neutrino 
are excluded. 



C. Atmospheric neutrinos 

The evidence for an anomaly in atmospheric neutrino experiments was pointed out by 
the Kamiokande Collaboration @,^] and 1MB Collaboration || using the water-Cherencov 
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experiments. More recently, SuperKamiokande Collaboration []19| , p0| reports the more pre- 
cise results on anomaly in atmospheric neutrino. We tabulate these results in Table [TI|. 
That the double ratio, R(ji/e) = i? expt (/z/e)/i?Mc(/V e )) * s l ess than 1 is the atmospheric 
neutrino's anomaly. In Table [HJ sub-GeV experiments detect the visible-energy less than 

I. 33GeV, and in the second column (total exposure), left number represents the detector 
exposure in which fully contained events are detected and right numbers partially contained 
events. 

The ratios -R exp t(/V e ) an d -RMc(/V e ) are defined as 

R t I \ _ / e n(Efj.) a n(.Eg, E^)F a {E a ,0 a )P{u a — » u^)dE^dE a d0 a 

expt W 6) ~ Eal e e (E e )a e (E Q ,E e )F a (E a ,0 a )P(u a - u e )dE e dE a d0 a ' [ ^ 

ft (u/ e ) — [ e ^{E il )a IJL {E m Efj)F a (E a , a )dE fl dE a d0 a (20b) 
Y, a f e e(E e )a e (Eg, E e )Fa(Ea,9a)dE e dE a d9g 

where the summation J2 a are taken in /i, e neutrino and /i, e untineutrino. €p(Ep) is the 
detection efficiency of the detector for /3-type charged lepton with energy Ep, ap is the dif- 
ferential cross section of up and F a (E a , 9 a ) is the incident u a flux with energy E a and zenith 
angle 9 a . P(u a — > up) is the transition probability Eq. (10) and it depends on the energy E a 
and the distance L which depends on zenith angle 9 a as L = JJr + h) 2 — r 2 sin 2 9 a — r cos9 a , 
where r is the radius of the Earth and h is the altitude of production point of atmospheric 
neutrino. 

Although informations of F a (E a ,9 a ) etc. are given in Refs. [p5|-|2"T|j, we use the MC 
predictions for f a (E a ,9 a ) = E« I t^E^a^Ea, E^)F a (E a , 9 a )dE fl in Ref. [§ for sub-GeV 
experiment and Ref. [0 for multi-GeV experiment. Explicit E a dependence of f a (E a ,9 a ) 
are shown in Appendix B. Since P(u a — > u^) and P(u a — > i/ e ) are the functions of 
(Am1 2 , Am% 3 ,6i2,9i3, 023, L,E), the double ratio R(/i/e) which is integrated in neutrino 
energy E and zenith angle (related to L) is the function of (Amf 2 , Am| 3 , ^12, ^13, ^23)- We 
estimate the R(fi/e) fixing the parameters (5m 2 2 , sm 2 20) on the allowed values Eq. (15) 
predicted by the solar neutrino experiments; 5m\ 2 = 3 x 10~ 5 eV 2 , sin 2 2^12 = 0.7 which cor- 
responds to large angle solution and 5m\ 2 = 10~ 5 eV 2 , sin 2 2# 12 = 0.005 which corresponds 
to small angle solution. 

In Fig. [|, we showed the contour plots of double ratio R(/i/e) on tan 2 13 — tan 2 #23 plane 
for various Am 23 . Contoure lines correspond the upper and lower values of R(fi/e) in Table 

II. We showed the plots of sub-GeV experiment in Fig. f|(a)-(d) and plots of multi-GeV 
one in Fig. [|(e)-(h), and in these figures solid lines denote the large angle solution plots 
and dotted lines the small angle solution plots. In Fig. (a)-(c), the allowed regions are 
surrounded by two solid lines (or dotted lines) and in Fig. (d), by two outer lines (or dotted 
lines) and inner solid line (or dotted line). In Fig. (e)-(h), dotted lines are close in solid 
lines. In Fig. (e)-(g), allowed regions are surrounded by two outer solid lines (or dotted 
lines) and inner solid line (or dotted line) and in Fig. (h), by two solid lines (or dotted 
lines) . 

Observing the allowed regions obtained terrestrial neutrino experiments Fig. ^| and the 
present allowed regions obtained atmospheric neutrino experiments, we obtain the allowed 
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regions satisfying all experiments. First, we show the allowed regions of (#13, 623) for the 
large angle solution and small angle solution in sub-GeV experiment: 



large angle solution 
for AmS^ = 20eV 2 no allowed region 

for Am| = 2eV 2 (0 13 < 4°, 18° < 23 < 26°) 

for Am 2 23 = 0.2eV 2 (0 13 < 4°, 18° < 23 < 62°) 

for Am 2 = 0.02eV 2 I 



for Am 2 3 = 0.002eV 2 



(013 < 12°, 17° < 23 < 38°) 
(013 < 12°, 47° < 23 < 63°) 

(013 < 29°, 27° < 23 < 52°) 



small angle solution 
no allowed region 

(013 < 4°, 24° < 23 < 26°) 
(013 < 2°, 64° < 023 < 66°) 
(0 13 < 4°, 24° < 23 < 66°) 
(013 < 14°, 23° < 023 < 41°) 
(013 < 14°, 49° < 023 < 67°) 

(013 < 34°, 34° < 23 < 56°) 



(39° < 0i3 < 59°, 76° < 23 ) [ (34° < 13 < 56°, 56° < 23 ) 
second, in multi-GeV experiment: 



for Am 2 3 = 20eV 



for Am| 3 = 2eV 



small angle solution 
no allowed region 

(0 13 < 4°, 24° < 23 < 26°) 
(013 < 2°, 64° < 23 < 66°) 
for Am 2 3 = 0.2eV 2 (0 i3 < 4°, 24° < 23 < 66°) (0 i3 < 4°, 24° < 23 < 66°) 
for Am 2 3 = 0.02eV 2 (0 13 < 12°, 41° < 23 < 49°) (0 13 < 12°, 41° < 23 < 49°) 
for Am 2 3 = 0.002eV 2 (22° < 13 < 68°, 45° < 23 ) (22° < 13 < 68°, 45° < 23 ) 



large angle solution 
no allowed region 

(013 < 4°, 24° < 023 < 26°: 

(0i 3 < 2°, 64° < 23 < 66°: 

9i 3 < 4°, 24° < 23 < 66°) 



(21a) 
(21b) 
(21c) 
(21d) 

(21e) 



(22a) 
(22b) 

(22c) 
(22d) 
(22e) 



If we combine the LSND experiment with above terrestrial data, allowed regions are 
restricted as follows, 
sub-GeV case: 



large angle solution 



for Am 2 3 = 2eV 2 (2° < i3 < 4°, 18° < 23 < 26°) < 



for Am 2 ^ = 0.2eV 2 no allowed region 



small angle solution 

(2° < 0i 3 < 4°, 

24° < 23 < 26°) 

(0.8° < 0i 3 < 2°, 

64° < 23 < 66°) 
no allowed region 



(23a) 



(23b) 



multi-GeV case: 
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for Am 2 23 = 2eV 2 



large angle solution 

(2° < 13 < 4°, 

24° < 23 < 26°) 

(0.8° < 9 13 < 2°, 

64° < 9 23 < 66°) 
for Am 2 3 = 0.2eV 2 no allowed region 



small angle solution 

(2° < 13 < 4°, 

24° < 23 < 26°) 
(0.8° < 0i3 < 2°, 

64° < 23 < 66°) 
no allowed region 



(24a) 



(24b) 



Although there are allowed regions satisfying both terrestrial and atmospheric exper- 
imental restrictions, all of these solutions do not satisfy the zenith angle dependence of 
R(fi/e) for atmospheric neutrino experiments. From resent SuperKamiokande experiments 
|H||2"0"|| , zenith angle dependence seems to be more definite than that obtained previously 
|||7|]: the double ratio R(fi/e) for sub-GeV experiment seems to decrease monotonically 
as zenith angle increases from = 0° to = 180° and R(fi/e) for multi-GeV experi- 
ments seems to decrease as zenith angle increases from = 0° to = 180°. Among the 
allowed solutions obtained above, Eqs. (21) and (22) (or (23) and (24)), the large angle 
(sin 2 20i2 ~ 0.7) solutions with 2 3 < 45° satisfy the monotonic decreasing of the R(fi/e) 
in sub-GeV neutrino experiment. This solution also decreases in multi-GeV experiment as 
zenith angle increases. 

In Fig. ||(a), we showed the zenith angle (cos0) dependence of R(fi/e) in sub-GeV exper- 
iment on the typical parameters for large angle solution (Am 2 2 = 3 x 10 _5 eV 2 , sin 2 20 i2 = 
0.7, Am 23 = 0.2eV 2 , 0i3 = 4°, 23 = 25°) by solid curve and small angle solution 
(Am 2 2 = 10" 5 eV 2 , sin 2 20 12 = 0.007, Am 23 = 0.2eV 2 , i3 = 4°, 23 = 25°) by dotted curve. 
Experimental data is referred to SuperKamiokande p0[ . Fig. |5|(b) represents the zenith an- 
gle (cos0) dependence of R(fi/e) in multi-GeV experiment on the parameters for large angle 

0.7, Am 23 = 0.2eV 2 , 13 = 4°, 23 = 30°) (solid 



10~ 5 eV 2 , sin 2 20i2 = 0.007, Am 23 = 0.2eV 2 , 
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solution (Am 2 2 = 3 x 10- 5 eV 2 , sin 2 29 12 
curve) and small angle solution (Am 2 2 
4°, 23 = 30°) (dotted curve). 

We summarize the solution satisfying the solar, terrestrial and atmospheric experiments 
as follows; 



for Am 2 3 = 2eV 2 



Am 2 l2 



4 x 10~ 6 - 7 x 10~ 5 eV 2 , 
13 < 4°, 24° < 23 < 26°, 



sin 2 20i2 = 0.6-0.9, 



(25a) 



for Am 2 3 = 0.2eV 2 

Am 2 12 = 4 x 10- 6 - 7 x 10- 5 eV 2 , sin 2 20 i2 = 0.6 - 0.9, 
013 < 4°, 24° < 23 < 45°. 

2 
23 



(25b) 



for Am 2 = 0.02eV 2 
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Am 2 12 = 4 x 1(T 6 - 7 x l(T 5 eV 2 , sin 2 29 l2 = 0.6 - 0.9, 



(25c) 



#i3<12°, 9, 



23 



~ 40°. 



If we include the LSND experiment in terrestrial ones, solutions, Eqs. (25a) and (25b), are 
favoured. 



We analyzed the solar, terrestrial and atmospheric neutrino experiments using 
the three-flavor neutrinos framework and got the allowed regions of the parameters 
(Am 2 2 , sin 2 26 l 1 2, Am 23 , 9 13 , 9 23 ). In solar neutrino experiments, we got the large an- 
gle solution (Am 2 2 , sin 2 2#i 2 ) = (4 x 10~ 6 - 7 x 10~ 5 eV 2 , 0.6 - 0.9) and small angle 
solution (3 x lO" 6 - 1.2 x 10" 5 eV 2 , 0.003 - 0.01) for 9 13 = 0° - 20°. When 9 13 in- 
creases from 25° to 50°, large angle and small angle solutions merge. From the terres- 
trial and atmospheric neutrino experiments (not considering the zenith angle dependence), 
we got the allowed regions: (#13 < 4°, 24° < 9 23 < 26°) for large angle solution and 
Am 2 3 = 2eV 2 , (9 13 < 4°, 24° < 9 23 < 26°) and (0 13 < 2°, 64° < 9 23 < 66°) for 
small angle solution and Am 23 = 2eV 2 , (9 13 < 4°, 24° < 9 23 < 62°) for large angle so- 
lution and Am^ 3 = 0.2eV 2 , (9 13 < 4°, 24° < 9 23 < 66°) for small angle solution and 
Am 2 23 = 0.2eV 2 , (9 13 < 12°, 9 23 ~ 40°, 48°) for large angle solution and Am 2 23 = 0.02eV 2 , 
(9 13 < 14°, # 23 ~ 41°, 49°) for small angle solution and Am 2 3 = 0.02eV 2 . When the 
zenith angle dependences are considered, the allowed regions are restricted to (Am 2 2 = 
4 x 10~ 6 - 7 x 10~ 5 eV 2 , sin 2 29 12 = 0.6 - 0.9, 9 13 < 4°, 24° < 9 23 < 26°) for Am 2 3 = 2eV 2 , 
(Am 2 2 = 4 x 10" 6 - 7 x 10~ 5 eV 2 , sin 2 29 12 = 0.6 - 0.9, 9 13 < 4°, 24° < 9 23 < 45°) for 
Am 2 3 = 0.2eV 2 and (Am 2 2 = 4 x 10" 6 - 7 x 10" 5 eV 2 , sin 2 29 12 = 0.6 - 0.9, 9 13 < 12°, 9 23 ~ 
40°) for Am| 3 = 0.02eV 2 . If we include LSND experiment in terrestrial experiments, allowed 
regions are restricted to (Am 2 2 = 4 x 10~ 6 - 7 x 10~ 5 eV 2 , sin 2 29 12 = 0.6 - 0.9, 2° < 9 13 < 
4°, 24° < 9 23 < 26°) and Am 2 3 = 2 - 0.2eV 2 . 

Finally, we present the neutrino mixing matrix Eq. (2) numerically for the allowed 
solutions Eqs. (25a), (25b) and (25c); 



IV. DISCUSSIONS 



for Am 2; 



2eV 2 



0.81 <-> 0.90 



0.43 <-> 0.58 



0.0 <-> 0.07 



-0.39 <-> -0.56 0.71 <-> 0.83 0.41 <-> 0.44 



(26a) 



for Am 2 = 0.2eV 2 



0.12 <-> 0.26 -0.33 <-> -0.42 0.90 <-> 0.91 



J 



11 



u 



for Am 2 3 



U 



0.81^ 0.90 0.43^0.58 0.0^0.07 

-0.30^-0.56 0.54^0.83 0.41 0.71 

0.12 <-> 0.41 -0.33 «-> -0.64 0.71 «-> 0.91 
= 0.02eV 2 

0.79 0.90 0.42 <-> 0.58 0.0 «-> 0.24 

-0.33 «-> -0.57 0.53 <-> 0.69 0.62 «-> 0.64 

0.11 «-> 0.38 -0.52 <-> -0.66 0.74 <-> 0.77 



(26b) 



(26c) 



APPENDIX A: NUMERICAL ANALYSIS FOR R IN SOLAR NEUTRINO 

The ratio R of the detected e neutrino flux to the expected e neutrino flux deduced from 
SSM is expressed as 



R 



Je::: 

where E is the neutrino energy and f(E) is the summation of products of fi(E) and cf 



(Al) 



f{E)=Y.cfh{E)- 



(A2) 



fi(E) represents the neutrino flux produced in the reaction of type i, where i = 
pp, 13 N, 15 0, 7 Be, 8 B and pep, at the center of the sun. cf is the detector sensitivity 
of detector D for the neutrino flux of type i. From the SSM [p3,2^|, neutrino fluxes at 1AU 
are expressed as the function of E (numerical number without dimension of neutrino energy 
measured in the unit MeV) as 

fpp(E) = £ L88 (4.358 - 0.0946E - 24.48E 2 ) x 10 12 , (E < 0.420MeV) (A3a) 

f Be i{E) = 8 x 10 8 , (0.384 - 0.25MeV < E < 0.384 + 0.25MeV) (A3b) 

f N (E) = £ L88 (4.194 - 4.541E + 0.870E 2 ) x 10 9 , (E < 1.199MeV) (A3c) 

f Be 2(E) = 4 x 10 9 , (0.862 - 0.5MeV < E < 0.862 + 0.5MeV) (A3d) 

fo(E) = £ L88 (12.74 - 7.919E + 0.328E 2 ) x 10 s , (E < 1.732MeV) (A3e) 

fpep(E) = 1.5 x 10 8 , (1.442 - 0.5MeV < E < 1.442 + 0.5MeV) (A3f) 

f B (E) = £ L88 (81211 - 11500£ - 407.1E 2 ). (E < 14.02MeV) (A3g) 

We estimate the detector sensitivity which is normalized as c^ = 1 using the Bahcall's result 
(Table | ||, 
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Ga experiment : c! 



Ga 
VP 



9.62 x 10~ 4 , c% a = 2.73 x 10~ 3 , c% a elBe2 = 3.61 x 1(T 3 , 



cg a = 4.30 x 10" 3 , c% 



8.93 x 1(T 3 , c% a 



1.0, 



CI experiment : c cl 



CI 



vv 



0, d N 



4.82 x 10-, c Bel 



r 1 



7.64 x 10- 4 , cg p 



Kamiokande experiment : c; 



Karri 
PP 



C 



Kara 
pep 



0, c B 



Kara 



= 0, c 
1.0. 



1.44 x 10~ 3 , c c B l 

Kara 



11 =5.08x10" 
1.0, 



Kara 
N 



0, c 



Bel 



0, C 



Kara 
Be2 



0, c, 



Kara 
O 



0. 



Numerical values used for parameters at the center of the sun are as follows: 

A = 2V2G F N e E = 2V2G F (Y e /m n )pE, 

G F = 1.17 x 10~ 23 eV~ 2 (Fermi constant), 

Y e = 1/2 (the number of electrons per nucleon), 

m n = 939 MeV (the nucleon mass), p = 156g/cm 3 (the density), 



(A4a) 
(A4b) 

i 

(A4c) 
(A5a) 



1 dN e 
W P ~dx 



2.91 x 10 



-15 



(A5b) 



APPENDIX B: EXPLICIT E a DEPENDENCE OF f a {E a ,9 a ) 

We take the E a dependence of f a (E a ,8 a ) from the results estimated in Monte Callro 
calculation. We neglect the 9 a dependence of f a (E a ,9 a ). For the sub-Gev experiments, we 



take it from Ref. [[J as 

f Vft (E, 6) = A S (-2.89E" 2 - 55 + 22.9E- 1 - 55 - 12.7£- Q " 55 ), (Bla) 

f Ue (E, 6) = A s (-1.18£- 2 - 55 + 9.98E- 1 - 55 - 3.90£T a55 ), (Bib) 

and for the multi-GeV experiments from Ref. ]7] as 

f Ufi (E, 6) = A m (42.1 - 383£" 2 + 374E" 1 - 2A2E + 0.0403E 2 (B2a) 

-0.000206E 3 ), fully-and partially-contained events (B2b) 

f Ue (E, 9) = A m (-9.54 - 225E- 2 + 2A9E' 1 + 0.156E - 0.000658E 2 (B2c) 

-2.42 x 10~ 6 £ 3 ), fully-contaiened events (B2d) 



where E is the numerical number of neutrino energy measured in unit GeV, and N s and N, 
are constants with the unit: number of events/GeV. 
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FIGURES 

FIG. 1. Contour plot of R for #13 = on sin 2 29\2 — Am 2 2 plane. Each two contours denoted 
as Ga, CI and Kam corresponds to the upper and lower value of R for Ga, CI and Kamiokande 
experiment, respectively. Two solutions are denoted as common areas enclosed by each two contours 
of Ga, CI and Kam. 

FIG. 2. The plots of allowed regions of the combined Ga, CI and Kam experiments using the 
X-square, where the solid thin, solid thick and dotted lines define the regions allowed at 99%, 95% 
and 90% C.L., respectively. Figs, (a) - (h) show the allowed regions for #13 = 0° — 50°. 

FIG. 3. The plots of allowed regions on #13 — #23 plane determined by P of terrestrial — ► u T , 
— > v e , — > v e and 9 e — > v e experiments. Allowed regions are corners, left and right hand 
sides surrouned by curves. Curves represent the boundary of 90 % C.L. of P. Ato 2 2 and sin 2 29\2 
are fixed as 10" 5 eV 2 and 0.8, respectively. Amjjg is fixed to 20eV 2 (Fig. 3(a)), 2.0eV 2 (Fig. 3(b)), 
0.2eV 2 (Fig. (c)) and 0.02eV 2 (Fig. (d)). Dotted lines show the allowed regions restricted by the 
LSND data. 

FIG. 4. The plots of allowed regions of R(fi/e) of atmospheric neutrinos for various values of 
Am 23 on tan 2 #13 — tan 2 #23 .Figs, (a)-(d) show the plots of sub-GeV experiments and Figs, (e)-(h) 
the ones of multi-GeV experiments. We fix the parameters (Am 2 2 , sin 2 29) to be (3 x 10 _5 eV 2 , 0.7) 
corresponding to the large angle solution (solid lines) and (10~ 5 eV 2 , 0.005) corresponding to the 
small angle solution (dotted lines). 

FIG. 5. The zenith angle dependence of R(fi/e) in atmospheric neutrino experiment. Ex- 
perimental data is referred to Ref. [20]. Fig. (a) represents the sub-GeV experimental case: 
solid curve is given on a typical parameters for large angle solution (Am 2 2 = 3 x 10 _5 eV 2 , 
sin 2 20i 2 = 0.7, Am 23 = 0.2eV 2 , 6>i 3 = 4°, # 23 = 25°) and dotted curve on small angle solu- 
tion (Am 2 2 = 10" 5 eV 2 , sin 2 29 12 = 0.007, Am 23 = 0.2eV 2 , 6 13 = 4°, <9 23 = 25°). Fig. (b) 
represents the multi-GeV experimental case: solid curve is given on parameters for large angle 
solution (6*i3 = 4°, 623 = 30°) and dotted curve on small angle solution (#13 = 4°, 623 = 30°). 
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TABLES 

TABLE I. Individual neutrino contribution (averaged on the cases with and without diffusion) 
calculated by Bahcall(1995) 



lie U. LI LLlKJ 0\J U.I be 






VP 
pep 
7 Be 
8 B 

13 N 
15Q 


0.00 
0.23 
1.17 
6.4 
0.09 
0.3 


70 

3 

35 
14 
3 
5 


total 


8.2 


132 


TABLE II. The results of atmospheric experiments predicting the anomaly. Sub-GeV experi- 
ments detect the visible-energy less than 1.33GeV, and in the second column (total exposure), left 
number represents the detector exposure in which fully contained events are detected and right 
numbers partially contained events. 


Experiments 


Total exposure(ktyr) 


i?(/i/e) (double ratio) 


Kamiokande(sub-Gev) || 
Kamiokande(multi-GeV) |7j 
1MB § 

SuperKamiokande(sub-GeV) [p0| 
SuperKamiokande(multi-GeV) |20|| 


4.92 
8.2, 6.0 

7.7 
20, 18 
20, 18 


0.60+o;°g ±0.05 
0.57+°;°® ±0.07 
0.54 ±0.02 ±0.07 
0.63 ±0.03 ±0.05 
0.60 ±0.06 ±0.07 
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